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Abstract 

If dark matter is made up of Weakly Interacting Massive Particles, the annihilation of these 
particles in halos results in energy being released, some of which is absorbed by gas, causing partial 
ionization and heating. Dark matter annihilation may result in partial ionization and gas heating at 
high redshifts, even before the formation of the first stars. It is shown that early ionization results 
in a transfer of power to higher multipoles in the large angle CMB polarization power spectra. 
Future CMB experiments may be able to place constraints on certain light dark matter models. 
We also investigate the effect of gas heating on the expected H21 cm power spectrum. Heating 
by particle annihilation results in a decrease in the amplitude of the H21 cm power spectrum as 
the gas temperature T becomes comparable to the CMB temperature T 7 , and then an increase 
as T > T 7 . The result is a minimum in the power spectrum at the redshift for which T ~ T~. 
Only certain models (low particle masses ~ 10 GeV, or favorable halo parameters) show this effect. 
Within these models, observations of the H21 cm power spectrum at multiple redshifts can help 
us obtain constraints on dark matter particle and halo properties. 
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I. INTRODUCTION. 



It is widely accepted on the basis of observations of the cosmic microwave background, 
gravitational lensing, galaxy clusters, rotation curves, etc, that ~ 80% of the matter in 
the Universe is in the form of a collisionless, non-baryonic component, commonly referred 
to as dark matter. Axions and Weakly Interacting Massive Particles (WIMPs) are two of 
the leading candidates for dark matter. Axion detection experiments include ADMX pQ, 
CAST [2], Germanium detectors [3], etc. WIMPs may be detected in a variety of ways, both 
directly and indirectly. Direct detection experiments include DAMA jJJ, CDMS [5], Zeplin 
[6], Edelweiss [7j, etc. Indirect detection experiments include various accelerator and gamma 
ray search experiments. WIMPs are also relevant to various astrophysical phenomena such 
as the formation of early stars [S], stellar evolution [H], etc. 

Supersymmetry provides a stable, neutral dark matter candidate called the neutralino. 
Neutralinos do not decay into standard model particles, but they annihilate in pairs, releasing 
photons, neutrinos, and charged particles. The energy released by particle annihilation 
is proportional to the square of the dark matter density and hence, most of the energy 
is produced by annihilation in dark matter halos. The earliest WIMPy halos form at a 
redshift z ~ 60, with masses M ~ 1O~ 6 M [ID], and are thought to be compact and 
abundant. Accretion and mergers result in more massive halos, in a bottom-up approach. 
We restrict ourselves to the case of thermally produced WIMPs with an annihilation cross 
section determined by the relic density today. 

Observations of the CMB polarization complement information obtained from the temper- 
ature anisotropy spectrum. The absence of strong Ly-a absorption lines (the Gunn-Peterson 
test [llj) in the spectra of quasars for z < 6 implies that the Universe is nearly fully ionized 
today. The WMAP measured optical depth may be used to place constraints on reionization 
models. Sudden reionization models with full ionization at redshifts below z* ~ 11 result 
in a peak in the CMB polarization power spectrum at low multipoles (large angles). The 
polarization power spectrum at large angles is sensitive to the assumed reionization model. 
Future experiments such as Planck [12] and CMBPol [13] may be used to distinguish between 
different reionization models [Hj. In this article, we compare the dark matter reionization 
scenario with the simple case of the sudden reionization model. 

Another powerful probe of the epoch of reionization, and the Universe prior to reioniza- 
tion, is the Hydrogen 21 cm line. Observations of the fluctuations in the H21 cm differential 
brightness temperature serve as a probe of the Universe at redshifts z < 200 [15J. Before 
the formation of the first luminous sources, the H21 cm signal is expected to be negligible 
at redshifts z ^ 30. This scenario is altered if dark matter annihilation provides a source of 
heating. Future experiments such as LOFAR [16J and the Square Kilometer Array [IT] aim 
to study the fluctuations in the H21 cm temperature, and provide important information 
regarding reionization and the high redshift Universe. Several authors have studied the im- 
pact of particle annihilation on reionization [TSJ [02 EH ED E21 E3] and the Hydrogen 21 cm 
signal [231 EU ESI EE]. In this article, we investigate possible signatures of dark matter an- 
nihilation in the CMB and Hydrogen 21 cm power spectra at redshifts before the formation 
of the first luminous objects, i.e. between z ~ 60 and z ~ 10. 

In Section II, we consider dark matter halos fitted with an NFW profile, and calculate 
the luminosity of the halos. We then compute the energy absorbed by gas at any redshift z 
due to particle annihilation in dark matter halos at redshifts z' > z. Some of the absorbed 
energy results in ionization and heating. We include the effect of gas heating in computing 
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the evolution of the ionized fraction, thus improving upon the calculation in [20] . 

In Section III, we calculate the optical depth due to scattering of free electrons with CMB 
photons. Early ionization results in more scattering of free electrons with CMB photons, and 
hence a larger polarization signal. However, the partial ionization at high redshifts requires 
a smaller ionization at low redshifts in order to keep the optical depth constant. Thus the 
large angle polarization power spectrum is modified. The effect is small, but detectable for 
small particle masses (m x ~ 2 GeV), and favorable halo parameters. 

In Section IV, we study the Hydrogen 21 cm differential brightness temperature. We 
show that the increase in gas temperature due to particle annihilation at high redshifts can 
have an observable effect on the power spectrum of fluctuations, for certain light dark matter 
models (m x ~ 10 GeV). As the gas temperature T increases and becomes comparable to the 
CMB temperature T 7 , the amplitude of fluctuations in the H21 cm differential brightness 
temperature decreases and reaches a minimum. The amplitude of fluctuations then increases 
as T > T 7 . We compare the power spectrum to the standard scenario that includes no dark 
matter heating. Finally, we present our conclusions. 



II. PARTICLE ANNIHILATION IN DARK MATTER HALOS. 

Consider dark matter halos fitted with a Navarro- Frenk- White (NFW) density profile |27j : 

P(r) = 2- (!) 

(r/r s ) [1 + r/r s ] 

p s and r s are constants. We define the concentration parameter C200 = ?"2oo/?V ^200 is 
the radius at which the enclosed mean density equals 200 times the cosmological average 
measured at the redshift zp(M) at which the halo formed: 

3M , 

= p{M) = 200p c fi m [1 + z F (M)f . (2) 
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M = M(r2oo) is the mass of the halo. We set the mass in dark matter to be Md m = 
(fidm/^m) M. p c is the critical density, and Q m is the matter fraction. C200 is thus measured 
at the time of halo formation. In this article, we treat C200 as a free parameter, independent 
of the halo mass. Particle annihilation in dark matter halos results in a luminosity (energy 
per unit time) per photon energy = dL/dE^ given by [28] : 

dL _ dN y ^ (a a v) f^ A _ 2n2fr) 



dE^ dE^ 



£Lx^4 i dr^r 2 p 2 
2 mi J 



dN 7 (a & v) M p /fi dm \ 



E 

dE y 7 2m? x 3 Va 



2 



) /(caoo). (3) 



m x is the particle mass. dN 1 /dE 1 = m" 1 dN^/dx is the number of photons released per 
annihilation per photon energy. Following |29j for mixed gaugino-higgsino models , we choose 
dN 1 /dx to have the form a e~ bx /x , where a and b are constants for a particular annihilation 
channel. Averaging over the channels considered in we find a = 0.9,6 = 9.56. (cr a t>) is 
the annihilation cross section of the WIMPs times the relative velocity, averaged over the 
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velocity distribution. We set (a a v) = 3 x 10 26 cm 3 /s, assumed independent of v (i.e. s-wave 
annihilation). For the NFW density profile Eq. , we have 



/(C200) 



4)0 [1-(1 + C200)- 3 ] 

3 [ln(l + c 20 o) - c 20 o(l + C200) 



(4) 



/(c 2 oo) takes values 18, 45, 150, and 610, for c 2 oo — 2.5, 5, 10, and 20 respectively. In this 
paper, we restrict ourselves to the density profile of Eq. [I]. The luminosity can be much 
larger for halo profiles that are more cuspy. For density profiles with p falling off faster than 
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, the luminosity is divergent unless a lower cutoff radius is imposed [20J. 



A. Interaction of photons with gas atoms. 



Consider dark matter annihilation at a location s'. Some of the photons released at s' 
travel to a location s where they scatter with gas atoms, resulting in ionization and heating. 
The probability of scattering between s and s + Ss is given by 5s cr(E y ) n [1 + z(s)} 3 where 



Ss = cSt = -cSz/ [(1 + z)H(z)} . 



(5) 



n(l + z) 3 is the physical number density of atoms at z. The fraction of photons that reach 
s having been released at s' is 



K(S, S', Ery) 



k(z, z', E y ) 



exp 



- exp 
2A 

'T 



a(E 7 )n [ S ds"[l + z(s")f 

J s' 

fa{Ej) [(1 + Z'f 2 - (1 + Zf' 2 



A is the dimensionless quantity 
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and f a (Ey) is the energy dependent correction to the Thomson cross section <tt (see for e.g., 



Uv) 



2(1 +y) | log(l+2y) 



2(1 + y) 2 2(l + 2/)log(l + 2|/) , 2 



1 + 2yf 



y 



l + 2y y 2 (l + 2y) 



y 3 
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where y = E^/m e and m e is the electron mass. For energies E y ^> 13.6 eV, we may ignore 
the bound structure of the atom, and consider scattering of photons with electrons. 

The comoving number density of dark matter halos at a redshift z is given by the Press- 
Schechter formula I3T1 



drib 
d~M 



[z,M) 



2 p c n m S c (l + z) da h 



vr M 



crt 



dM 



exp 



2al 



(9) 



o\ is the variance of halo fluctuations (see for e.g., [32]): 

dk k 3 P{k) 



W(kR) 



2tt 2 



\W(kR)[ 



(kRY 



[sin kR — (kR) cos kR] . 



(10) 
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We use the code of Eisenstein and Hu [33] to compute the matter power spectrum P{k). 
P(k) is normalized so that a 8 = a^(R = 8/i _1 Mpc) = 0.8. S c is set to 1.28 |34j. Assuming 
density fluctuations grow as (1 + z)' 1 , the matter power spectrum at redshift z is P(z, k) = 
P{k)/{l + z) 2 . 
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FIG. 1: (a) shows the average absorbed energy per atom per unit time (E/n(l + z) 3 ) plotted as 
a function of redshift, for particle masses m x = 10,50, and 100 GeV, with C200 = 10. (b) shows 
E/E*. Peak absorption ~ 20% for m x = 10 GeV, and ~ 6% for m x = 100 GeV. 

The energy absorbed by gas at a redshift z is E(z, M min , m x , c 2 oo) given by (see [20] for 
details) 

E = H5^> /{<*>) f^)\ 1 + zf f _ dz , (1 + jylfi 
/ ( c 200 — lUJ V I'm / 



<>X J v^200 w y v "m ' jz f 

x f 1 dx af ^f a (m x x)K(z 1 z',m x x) f°° dM ^(z' , M) p(M) M. (11) 

J y I «/ M m ; n UikZ 



£p — Zp ( Mmi n) is that redshift at which cTh(M min , z-p) = 1 provided zp > z [10J. If not, we 
set z-p = z. We use the form given by [10J to estimate zp(M). The energy absorbed by the 
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FIG. 2: Different choices of the minimum halo mass M m ; n . E is small at late times because of the 
low halo and gas densities, and at very early times because halos are rare. The formation redshift 
is larger for smaller M mm . The particle mass m x = 50 GeV. 

gas contributes to ionization as well as heating. Monte Carlo simulations [35] suggest that 
about a third of the absorbed energy goes into ionization, while a third goes into heating. 
We therefore set the efficiency factors i] ion = ?7 heat = 0.3. Fig. [ija) shows the average energy 
absorbed per atom per unit time (E/n(l + z) 3 ), for particle masses m x = 10,50, and 100 
GeV. 

It is instructive to compare E to the simple approximation E*: 

E m (z) = J dM(l+zf^(z,M) J dE^iE^M), (12) 

which neglects the effect of propagation of photons, and unlike E, is a strictly local function 
of z. Fig. [jjb) compares E and E*. We see that for light dark matter particles (m x ~ 10 
GeV), about 20% of the energy released is absorbed at high redshifts, falling to pa 10% at 
lower redshifts. This is in agreement with the somewhat larger absorption fraction obtained 
by [IH] at z — 50, for very light dark matter with mass m x = 10 MeV. Assuming the 
efficiency factors rj lon = r^eat = 0.3, we find that ^5 0.3 x 0.2 = 6% of the generated energy 
goes into ionization, and an equal amount into heating. The maximum absorbed fraction 
is lower (pa 6%) for 100 GeV dark matter particles, with a useful fraction < 2% each for 
ionization and heating. Fig. [2] shows the effect of varying the assumed minimum halo mass 
M min . The minimum halo mass is set by the free streaming scale which depends on the 
particle physics of the model. For WIMP dark matter, we expect M min pa 1O~ 6 M [TU]. We 
set M m i n to 1O _6 M in the remainder of this article. We now calculate the evolution of the 
ionized fraction and the gas temperature. 
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B. Ionization of gas. 



The number of ionizations per unit volume per unit time, at a redshift z is given by 

I(z) = n (1 - x ion ) T] ioQ E(z). (13) 
fi is the inverse ionization potential 

'0.76 1 0.06 1 



0.07 eV" 1 . 



(14) 



.0.82 13.6 eV 0.82 24.6 eVJ 

where we assumed 76% H and 24% He by mass. Only singly ionized Helium is considered in 
Eq. (14), and x- mn is the ionized fraction. The number of recombinations per unit volume, 

(15) 



per unit time at z is given by 

R(z)=n 2 x 2 OQ (z)(l + z) 6 a(T). 
T is the gas temperature and a is the rate coefficient 

a(T) 

with an and au e given by [36J 

a H « 3.746 x 10~ 13 (T/eV) 
a He ~ 3.925 x 10- 13 (T/eV)-°- 6353 

x ion is obtained by solving the equation 



0.76 /rT1 , 0.06 

^«H(T) + — «He(T) 



(16) 



°- 724 cm 3 s- 1 



3 — 1 

cm s . 



I(z) - R{z) = -nH Jn m (l + z) 



jll/2 



dz 



where we used 



dz 



dtH Jn m (l + z 



.5/2 



l + -*(l+z)- 3 

1 6rri 



,1/2 



(17) 



(18) 



(19) 



and neglected the effect of dark energy for z > 10. In order to solve Eq. (18), we need to 
obtain an expression for T(z). 



C. Heating. 

In the absence of heating, the gas temperature falls off oc (1 + z) 2 . Dark matter annihi- 
lation results in a heating term (dT / dz\ eat 

3 , tt /7^~fi , nii/2 I dT\ 



7] hcat E{z) = --nk h H jQ m {l + Z) 11 ' 1 — , (20) 



2 v ' \dz) 



heat 



where is Boltzmann's constant. In the presence of free electrons, energy is transferred 
between the photons and the gas by Thomson scattering, resulting in a coupling term pU [37] 

Tj[z) — Tjz) x ion (21) 

t c {z) l + X ion + / H e' 
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where a is the radiation constant, and /he is the Helium fraction. t c (z) = 3m e c/9>aa^T*(z) rs 
1.44 Myr [30/(1 + z)] A . This coupling increases the gas temperature when T is smaller 
than the CMB temperature T 7 and decreases it when T > T 7 . The evolution of the gas 
temperature thus follows the equation 



dT 

dz 



2T 



1 + Z H y/TU(l + Zfl 2 



+ 



r 7 -r 



2?7heat£' 

3A; b rz(l + 2;) 3 1 + x ion + / He t c 



(22) 



Solving Eq. (22) together with Eq. (18) gives T(z) and x 
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FIG. 3: Evolution of the ionized fraction x lon with redshift, for different values of the concentration 
parameter 0200- (a) shows the case for m x = 10 GeV while (b) is plotted for m x = 50 GeV. The 
residual ionized fraction at a redshift z = 55 was chosen to be 10~ 4 . 

Figs. [3] and [4] show the evolution of the ionized fraction x\ OIi and the gas temperature T 
with redshift. We note from Fig. [3]that only models with small particle masses m x ~ 10 GeV 
and concentration parameters C200 ~ 10 are important from the point of view of reionization. 
The solid red line in Fig. [4] shows the gas temperature ~ (l + z) 2 in the absence of dark 
matter heating. The solid black line shows the CMB temperature ~ (l + z). The dashed 
lines show the gas temperature for different models with dark matter heating. Our results 
are in agreement with previous work [23]. In Fig. g}a), the gas temperature eventually 
becomes larger than the CMB temperature for all the models considered. In Fig. |4(b), T 
exceeds T 7 only for the models with C200 — 10 and 20. 



III. EARLY REIONIZATION. 

Let us now consider different dark matter models and study the epoch of reionization for 
a fixed optical depth. We then compare the expected CMB polarization power spectra for 
the different models. 
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FIG. 4: Evolution of the gas temperature with redshift. (a) shows the case for m x = 10 GeV, 
while (b) is plotted for m x = 50 GeV. The solid red line shows the gas temperature ~ (1 + z) 2 in 
the absence of dark matter heating. The solid black line shows the CMB temperature ~ (1 + z). 
The dashed lines show the evolution of the gas temperature for various concentration parameters. 



A. Optical depth and reionization redshift. 

The optical depth due to scattering of CMB photons with free electrons is given by 

T ( z ) = J Q cdz' (^-^j a T n e (z'). (23) 

Studies of quasar spectra have confirmed that the Universe is nearly completely ionized 
up to z ~ 6. Assuming H is ionized at z = 6, He is singly ionized at z = 6, and doubly ionized 
at z = 3 [38], we find r(z = 6) = 0.04. As this value if less than the WMAP measured [39] 
value of r m 0.087, we conclude that the Universe is at least partially ionized at redshifts 
z > 6. Let z* be the redshift below which x lon = 1. The simplest model of reionization is one 
in which Xi on (z) = for z > z*. In this model, the Universe is assumed to be instantaneously 
reionized at z = z*. Let us refer to this model as the "sudden reionization" scenario. While 
only an idealization, this model gives us a rough idea of the epoch of reionization. Then, we 



can rewrite Eq. (23) as 




r = 0.04+ / cdz — — a T n(z) + cdz — — a T n(z)xi on (z). (24) 



The sudden reionization scenario gives z* = 10.66 assuming only ionized Hydrogen beyond 
z = 6. When ionization by dark matter is taken into account, the value of z* is decreased as 
a^ion 7^ for z > z*. Fig. 5 shows 2 different reionization histories, giving the same value of 
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FIG. 5: Two different reionization histories. The solid (red) curve shows a sudden reionization 
scenario in which the Universe remains neutral until a certain redshift z = z*, below which x; on = 1- 
The blue curve takes into account ionization by dark matter, for the model with m x = 2 GeV and 
C200 = 10. z* is different from the previous case and xi on ^ for z > z* . Both curves give the same 
optical depth r = 0.087. 
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(b) m x = 2 GeV 
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TABLE I: (a) shows how depends on the assumed particle mass m x , with C200 set equal to 10. 
(b) shows the variation of the reionization redshift with the assumed concentration parameter 
C200 with the dark matter particle mass m x set equal to 2 GeV. In each case, the total optical depth 
is 0.087. The cases C200 = and m x = 00 correspond to no ionization by dark matter annihilation. 

t = 0.087. The solid (red) curve shows the sudden reionization scenario in which x ion = 
for z > = 10.66 and x ion — 1 for z < z*. The dashed (blue) curve shows an alternate 
model which includes dark matter ionization for z > z* = 6.85. x\ on = 1 for z < z*. m x = 2 
GeV and C200 = 10. 

Table 1 shows z* for different values of m x and 0200- is considerably lowered only 
for very light dark matter particles (~ 1 GeV), and for large concentration parameters 
(c 2 oo ~ 10). 
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NoDM ionization 
With DM ionization 




FIG. 6: CMB polarization power spectra, (a) shows the EE power spectrum, while (b) shows 
the TE power spectrum. The solid (red) curve is plotted for the sudden reionization scenario with 
no ionization by dark matter (the solid, red curve in Fig. 5). The dashed (blue) curve includes 
ionization by dark matter and corresponds to the dashed, blue curve of Fig. 5 (m x = 2 GeV, 
C200 = 10). 

B. Effect on CMB polarization spectra. 

The partial ionization of the Universe at high redshifts increases the scattering of CMB 
photons with free electrons, thus influencing CMB polarization. The increased scattering at 
high redshifts results in a transfer of power from low multipoles to higher multipoles in the 
large angle EE and TE polarization spectra, compared to the sudden reionization scenario. 
Fig. [6] shows the EE and TE polarization power spectra computed using the CAMB software 
[10], for the two reionization histories considered in Fig. [5j The solid (red) line corresponds 
to the sudden reionization scenario in which x ion = 1 for z < z* and for z > z*. The 
dashed (blue) line corresponds to the reionization history with DM annihilation taken into 
account, for a model with m x = 2 GeV, C200 = 10. x lori = 1 for z < z*. (compare with 
Fig. [5]). The scattering of electrons at high redshifts results in a reduced amplitude for 
the first peak at / ~ 4 and greater power at slightly higher I. For I > 30, there is very 
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little difference between the spectra. It may be possible to distinguish these curves with 
future experiments such as Planck and CMBPol [131 EH]- However, differentiating between 
the dark matter models and a model of gradual reionization by astrophysical sources (as 
opposed to the sudden reionization model) will be considerably more challenging. It was 
recently shown [22J that in certain dark matter models with a large value of (cr a v), ionization 
by dark matter annihilation is substantially larger. This would result in large modifications 
to the CMB polarization power spectrum as well. 



IV. THE HYDROGEN 21 CM LINE. 



Hydrogen 21 cm cosmology is an important tool to study reionization and the early Uni- 
verse [15]. The 21 cm differential brightness temperature (brightness temperature relative 
to the CMB) is given by (see for e.g., ET 



T b (z) » 27 mK 



'1 + 2 

10 



(1 - Xj, 



n 
n 



T 



H(z)/(l + z) 
dvu/dru 



no is the spatial average of n. T s is the spin temperature given by 



t: 
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1 + £c + £a 

Cc{z, Xion, T) is called the collisional coupling coefficient [25] : 

n(l + z) 3 r„ 



Aio T 7>0 (1 + z) 



X,\, 



1 



(25) 



(26) 



(27) 



Aiq is the Einstein coefficient for spontaneous emission. T* = E w /kh, with E w being the 
energy difference between the singlet and triplet levels. K e (T) and k h (T) are rate coefficients 
for collisions with electrons and Hydrogen atoms respectively. At low temperatures, the 
ionization fraction is very low, and we expect k h to be the dominant term. n e becomes 
relevant for larger values of x lon which are associated with higher T. We may neglect 
collisions with protons for large T [41J. dv\\/dr\\ is the gradient of the proper velocity along 
the line of sight. £ a is the Wouthuysen-Field (Ly-a) coupling term [25] : 



16n 2 cf ( 
27A 




dvJ{v)®(v) 



(28) 



Fig. [7] shows £ c and £ a for m x = 10 and 50 GeV, with C200 = 10. For redshifts z < 25, 
both coupling coefficients are important for light dark matter particles m x ~ 10 GeV. For 
heavier dark matter particles m x ~ 50 GeV, £ c 3> £ a . f a = 0.4162 is the oscillator strength 
of the transition. J{v) is the number of Ly-a photons per unit area, per time, per frequency, 
and per solid angle. Following [24], we make the assumption that one half of the energy 
that goes into collisional excitations results in the generation of Ly-a photons. With the 
simple ansatz that a third of the total energy absorbed goes into collisional excitations, 
we have a fraction rj a = 1/6 of the total energy resulting in Ly-a photons. $(^) is the 
line profile satisfying / dv$>(v) = 1. We will assume that $(^) is sharply peaked about 
v = u a ~ 2.5 x 10 15 Hz. J may be expressed as [24] : 



J 



VaC 

47T Hn 



E 1 

z) 3 / 2 hv a v a 



(29) 



12 




FIG. 7: Collisional (£ c ) and Ly-a (£ a ) coupling terms for particle masses m x 
C200 = 10. 



10 and 50 GeV. 



E(z, m x , c 20 o) is given by Eq. 11 Using Eq. 28 with f2 m = 0.258, h = 0.71, we find 



£ Q « 0.012 



21 \ 5 / 2 

1 + 2 



10- 20 eVcnrV 1 



(30) 



Using Eq. (26), we can rewrite Eq. (25) as 



T b = 27 mK 



'1 + z 
10 



T7 



n i + e 



#(*)/(! + *) 
dv\\/dr\\ 



(31) 



where £ = £ c + £ a . 

Fig. [8] shows the spatially averaged value of Tb for different values of the concentration 
parameter c 2 oo, f° r masses (a) m x = 10 and (b) 50 GeV. The solid (red) curve shows T b in 
the absence of dark matter heating (C200 — 0). Let us refer to this temperature as Tb(no 
DM). The gas temperature with no dark matter heating is labeled T(no DM). For some 
values of C200, the differential brightness temperature starts out smaller than Tj-,(no DM), 
but soon crosses this curve. This is caused by the gas temperature T being smaller than 
the CMB temperature T 7 , yet larger than T(no DM). Eventually T becomes larger than T 7 
(compare with Fig. [3J. For (C200 = 20, m x = 10 GeV), Tb is larger than Xb(no DM) for all 
z. Conversely, for (C200 = 5,m x = 50 GeV), Tb is smaller than Tb(no DM) for all z. We will 
explore the redshift variation of T b further when we discuss the power spectrum. 
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C200 = 
C200 = 5 

c 2 oo = 10 
c 20 o = 20 




20 30 40 

(b) m x = 50 GeV 



FIG. 8: H21 cm differential brightness temperature The solid (red) line shows Tb with no dark 
matter heating, while the dashed lines are drawn for models with C200 = 5, 10, and 20. Shown are 
the plots for (a)m x = 10 GeV and (b)m x = 50 GeV. 



A. Perturbations in the brightness temperature. 

Let us now consider perturbations in the differential brightness temperature Tb. Eq. 



(31) is sensitive to fluctuations in the baryon density, the ionized fraction, and the gas 
temperature. Let us perturb Xi on , n, and T, and expand to linear order: 



1 - Xi, 



1 >^ion,0 3'ion,0^x 



n = n + n 5 n 
T T T 

T T T T 



£ = £0 + 



<9£o x d£ <9£ / 
7-" n oc>n + t^t T d T + x ionfi d x 



(32) 



where Xi on fi,n , an d T are spatially averaged values, and S x ,8 n , and 5^ are fractional per- 
turbations in the ionized fraction, baryon density, and gas temperature respectively. £ is 
evaluated at (x ioni0 , n , T ). 

The fractional perturbation in the brightness temperature in Fourier space is then found 
to be 



8 21 (z,k) = 6 n (z,k) 



+ 6 T {z, k)p T {z) + 6 x (z, k)fi x (z), 



p n (z) + (h-k)' 

where n is the direction of the line of sight, k is the Fourier wave vector, and 

1 <91n£o' 



(33) 



1 



1 + £0 d In n 
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Ac 





1 - T^/Tq 

>^ion,0 



1 d\n£o 
l + £o d\nT 
1 «91n£ 



+ 



1 — ^ion,0 1 + £o 9 In 



(34) 



To compute the derivates of £ = £ c + £ a , let us first consider perturbations in the density of 
dark matter halos. The energy absorbed by the gas at any redshift z depends on the halo 



density for redshifts z' > z, as seen in Eq. (11). The integration over redshifts in Eq. (11) 



results in the perturbations averaging out as statistically, overdense regions are as likely to 
exist as underdense regions. As a first approximation, we may ignore perturbations in halo 
density. The fractional perturbation caused by the dark matter halos oc 1/ VN where iV is 
the number of halos in a mean free path volume. Hence this assumption is valid provided 
VN > 1. The mean free path of a high energy photon at redshift z ~ [n(l + z) 3 a] w 100 
Mpc for z = 25, where n is the comoving baryon number density, and we assumed a = <tt 
(although for high energy photons, a < o~t, which makes the mean free path larger). The 
number density of halos at redshift z is 



"haloUJ 



roc 

(1 + ^) 3 / 



dM 



drib 
dM 



(35) 



10 18 Mpc 



-3 



(10 12 Mpc 

13 



-3\ 



for M, 

-3 



mm 
13 



10 M (1M ), at z = 25. The corresponding mass 
density « 8 x 1O 13 M Mpc~ 3 (10 13 M Mpc -3 ), about 15% (2%) of the dark matter density 
at z = 25. dn^/dM is given by Eq. |9]. Thus y/N ~ 10 9 — 10 12 depending on M min . Using 



Eq. (27) and £ Q oc E oc n, we have: 



<91n£ 
9 Inn 
<91n£ 
<91nT 
<91n£ 

9 In X inr , 



x ion dn c /dT + (1 - Xio^dn^/dT 



x 



ion 

c 



K e + (1 - Xi, 



,K e + 1 



H ' 



(36) 



We show that for redshifts 2 < 30, /3 n is often the dominant term. 



B. The multi-frequency angular power spectrum. 

Let us expand S21 as a sum over spherical harmonics: 

821(2, h) = 53 a im( 2 ) ^m(n) 



(37) 



where 

a«m(2:) = 
and r is the comoving distance 



d 3 k 



r(z) = t[^(-2)] — I dz 



exp—ikr(k-h) 821(2, k) 



2c 



H(z') H y/Q T 



i-a + z) 



-1/2 



(38) 



(39) 
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where we ignored Qa- v(z) is the redshifted 21 cm line frequency 



u(z) 



1 + z' 



(40) 



where u = c/21.1cm = 1.42 GHz. Setting x = kr, and using Eq. (33) and the identities 



(k-hf = -4tt(— i) 



dx 2 



we find the following expression for ai m : 

d 3 k 



a-imty) = 4vr(-i) 



(2tt 



d 2 ji(x) 
dx 2 



+ 5 t Pt + 5 x p x ji(x) 



(41) 



(42) 



We can now construct the variance Ci defined as 



Q(u, Av)5u>6 mm i = (a lm {v) a* Vm ,{p')) = 16tt 2 (— 



d 3 k d 3 k' 



2tt 3 2tt 3 



Y* m {k)Y Vm ,{k') V{z,z',k,k') 

(43) 



V is given by 

V = {5 n {z)5l{z'))[(3 n {z)j l {x)-d 2 j l {x)/dx 2 ] [f3 n (z')j v (x') - d 2 j v {x')/dx' 2 

+ j t (x) [Uz') 3l ,(x')-d 2 3l ,{x')/dx' 2 ] [^ T (z)(5 T (z)5:(z / ))+/3 x (^)(5 x (^^(z / ))] 

+ Ji'(x') [PMMx) -d 2 Jl {x)/dx 2 ] [p T (z')(5 n (z)5* T (z')) +^')(^WAz'))] 
+ Ji{x)jv{x') [p T {z)fc{z')(8 T {z)8* T {z')) +/3 x (z)/3 x (z')(5 x (z)^(/)) 
+ fr(*)&CO<M*)£CO> + /3x(^)^t(^)(^(^)^(^))] . 

x' = fcV. Let us now make the approximation that \u' — v\ = Av <C z/ : 

c 



|r' — r| = Ar « y/l + z 



Av 



(44) 



(45) 



For small values of \Av/vq\, we may make the following approximations: 

(5i{z) « A(^) 

(fcC*,*)^*',*')) ~ (5i(-z,fc)^(^fc')> 

= {2n) 3 5 3 {k - k')Pij{z, k) = {2n) 3 5 3 {k - k^P^z, k) 



(46) 



where i and j may stand for the baryon number density, ionized fraction, or gas temperature 
fluctuations. Using the above approximations, setting x = kr,x' = kr', and the identity 



J dVLYi m {h) Y^ m ,{n) = 5i V 5 n 



(47) 
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we may simplify Eq. ( 43 ) : 

2 r°° 

Q(u,Av) = - dkk 2 x 

7T JO 



P im (z, k) {(3 n {z) 3l {x) - d 2 3l {x)/dx 2 } {MzMJ) - d 2 Jl {x ! )/dx 12 } 
+ {(3 T (z)P nT (z, k) + (3 x (z)P nx (z, k)} [jfa) {Pn(z)ji(x') - d%(x')/dx' 2 } 

+ j l (x') {p n (z) 3l (x)-d 2 3l (x)/dx 2 } 
+ ji(x)ji(x') {/^P T tO, k) + 2(3 T f3 x P Tx (z, k) + f3 2 x P xx (z, k)} 



(48) 



C. Fluctuations in x; nr , and T. 



To compute Ci, we first need to compute the various power spectra in Eq. (48). Since 



we are ignoring perturbations in halo density, the fractional perturbation in E is then 5e = 
8 n . Let us express the fractional perturbations in the ionized fraction (5 X ) and the gas 
temperature (8t) at redshift z as: 

5 x (z,k) = S x (z)5 n (z,k) 

S T (z,k) = S T (z)5 n (z,k), (49) 
i.e. we make the approximation that the k dependent part of 5 X and St follows the form 



of 5 n . This is motivated by the dependence of Eq. (11), on the baryon density. S x and 



St give us the time evolution of the perturbations. With this approximation, we consider 



perturbations in Eq. (18) and Eq. (22), to obtain the linearized equations: 
nx ion ay/l + z 



dS x 
dz 

dS T 



1 



1 

l + z 



,L\, 



1 dx 1DI 
dz 



+ 



E 



S x 1 



T 
T 



Si 



+ S^ 



dz HoV^K(l + fH C ) t c (l + Z )^ 

where we used x- lon <C 1. The power spectra then become 

PTT(z,k) = S 2 ( Z )P nn ( Zj k) 

PU^k) = S 2 (z)P nn (z,k) 

P nT (z,k) = S T (z)P nn (z,k) 

P nx (z,k) = S x {z)P nn (z,k) 

PTx(z,k) = ST(z)S x (z)P nn (z,k). 



1 dT 



l + z T dz 



(50) 



(51) 



Eq. (48) can be simplified to the form 
2 



Ci(u,Au) 



I roc 

, / dkk 2 P(k) X 

7T (1 + Z f Jo V ' 



(52) 



where (3{z) = 
P(k)/(l + z) 2 . 



f3 2 {z)j l {x)j l {x')-^z)\j l {x 

Pn(z) 



d 2 ji{x') . , , s &ji{x) 
+ 3i{x 



dx' 2 dx 2 J 

St{z)Pt(z) + S x (z)/3 x (z), and we assumed P l 



+ 



d 2 ji(x) d 2 j l (x') 



dx 2 
z,k) = 



dx' 2 
P{z,k) 
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FIG. 9: /3Tb ,o and /3 n Tb,o as a function of redshift z. The curves are nearly identical at low redshifts 
which implies that our results are not very sensitive to the form of 5t and <5 X . 



I = 500 




io~ 6 io~ 5 io~ 4 10- 3 



FIG. 10: Variation of Ci(v, Av) with Av/vq, at z = 25 (v(z) ~ 55 MHz). The decrease in power 
with Av may be used to separate the Hydrogen 21 cm signal from the large foregrounds. 



Fig. [9] shows /3Tb 5 o and fl n Tb,o for particle masses m x = 10 and 50 GeV, with initial 
conditions St = S x = 1 at z = 55. T^o is the spatially averaged value of Tb. The difference 
between the two curves is a measure of the importance of fluctuations in the ionized fraction 
and temperature. In the redshift range 10 < z < 30, the two curves are nearly identical, and 
thus our results are not very sensitive to fluctuations in the ionized fraction and temperature 
in this range, except when T«iT 7 . 



Fig. 10 shows the normalized multifrequency power spectrum Ci(u, Au) as a function of 



Av I v$. The amplitude of the power spectrum decreases quickly as Av is increased from zero 
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because the Bessel functions oscillate out of phase. For I = 500, the amplitude falls to 1% of 
the peak value when Av ~ 1.3 MHz. For multipoles I = 1000 and I = 5000, the amplitude 
falls to 1% of the peak value for Av »s 0.8 and 0.2 MHz respectively. For comparison, the 
spectral resolution of the LOFAR experiment ~ 10 kHz or better [42J. The strong frequency 
dependence of the multi-frequency angular power spectrum may be used to separate the 
signal from the much larger foregrounds which are expected to be highly correlated over 
such small frequency ranges [43] . 



D. Variation with redshift. 



Let us now consider the variation of the amplitude of the power spectrum with redshift 
or equivalently, with frequency v = uq/(1 + z). We consider here, redshifts that are well sep- 
arated from each other. We assume that the signal has been separated from the background 
and plot the peak value Ci(v,Av = 0) for each redshift z. 



(M 

c 

+ 




(a) m x = 10 GeV 
C200 = 5 



10 2 10 3 10 4 
I 

(b) m x = 50 GeV 
C200 = 5 



(c) No DM heating 



FIG. 11: Power spectrum C[(v, Av = 0) for different z. (a) and (b) show plots for m x = 10 GeV 
and 50 GeV respectively, with C200 = 5. (c) is plotted for the case of no heating by dark matter. 
In (a), the power spectrum decreases, and then increases, as we move from redshift z = 15 to 30, 
in contrast to (b) and (c). 



Fig. 11 shows the H21 cm power spectrum Ci(u, Av = 0) for redshifts z = 15,20,25, 



and 30. Plots (a) and (b) are shown for m x = 10 and 50 GeV respectively, for an assumed 
concentration parameter C200 = 5. (c) shows the case of no dark matter heating. In (a) the 
amplitude of the power spectrum of fluctuations decreases as z is increased from 15 to 25, 
and then increases from 25 to 30. The variation of C/ with z is due to the heating of the 
gas by dark matter annihilation. As the gas temperature increases, T approaches T 7 , and 
hence the differential brightness temperature decreases, resulting in a loss of power. The 
power spectrum reaches a minimum when T = T 7 and then increases again as T > T 7 (the 
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21 cm background now appears in emission). This behavior is in contrast with (c) which 
shows a steady increase in power as the redshift is varied from 15 to 30. Thus observations 
of 21 cm fluctuations at different redshifts may be used to identify light dark matter models. 
The curves in (b) do not show a minimum with change in redshift. This is because for 
C200 — 5 and m x = 50 GeV, the gas temperature T never rises above the CMB temperature 
(see Fig. 4). Hence this effect is only visible with light dark matter models, or with large 
concentration parameters. 



Fig. 12 shows the effect of variation with redshift, for I = 1000, and for C200 = 2.5, 5, and 



7.5. Also shown is the case C200 = which does not take dark matter heating into account, 
(a) is plotted for m x = 10 GeV, while (b) shows m x = 50 GeV. In (a), the power spectrum 
shows minima for C200 = 2.5,5, and 7.5. In contrast, in plot (b), we see a minimum (at 
z = 16.4) only with c 2 oo = 7.5. Thus, for model (b) we would expect the power spectrum 
to show a minimum at this redshift when m x = 50 GeV, C200 = 7.5. No minimum exists for 
C200 = 2.5 or 5. 




FIG. 12: The power spectrum Ci(u, Av = 0) for I = 1000 as a function of z, for C200 = 2.5,5, 
and 7.5. v(z) = vq/{1 + z). Also shown is the case C200 = which corresponds to no dark matter 
heating, (a) and (b) are plotted for m x = 10 GeV and 50 GeV respectively. 



V. INFLUENCE OF ASTROPHYSICAL SOURCES. 

So far, we have ignored the existence of stars and other luminous astrophysical sources, 
and considered heating solely due to dark matter annihilation. For redshifts z < 15, this 
assumption is bound to break down as luminous astrophysical objects heat up and ionize the 
Universe. However, we expect that dark matter annihilation remains the dominant heating 
source at redshifts z > 20. Let us now consider the effect of astrophysical sources on the 
brightness temperature. 
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A. Region around early stars 



Let us first consider the formation of the earliest stars. These stars are massive, hot, 
metal-free, and form in dark matter halos with masses Mh a i > 7 x 1O 5 M |44J. Each such 
star is surrounded by a hot, ionized region (radius i?;) which in turn is surrounded by a 
warm neutral region (radius R n ). After the star has exhausted its fuel in a short time 
Az « 0.1 (Ai/Myr) at z = 25, the gas surrounding the star cools adiabatically, while the 
ionized region becomes partially neutral. In this simple estimate, we ignore the formation 
of supernovae, black holes, etc. 

Let us consider stars with mass M sta r ~ 100 M and lifetime ~ few Myr. The size and 
temperature of the ionized and neutral regions surrounding Pop. Ill stars have been estimated 
in [43] . Following [45] , we set a threshold temperature of 10 4 K for the central ionized region 
up to Ri = 3 kpc, at an assumed formation redshift z s = 25. The gas temperature drops 
quickly to the cosmological average of ~ 10 K at R Q = 8 kpc. For simplicity, we approximate 
the temperature dependence on radius of [45] by a power law T(r, z = 25) ~ 10 4 K (Ri/ r) 7 - 13 . 
The differential brightness temperature averaged over the sphere surrounding a single star 
is 

(T b ) =M J ' dQ J dsT h {s,6), (53) 

where the integrals are over the solid angle subtended by a single sphere when observed by an 
instrument, and the line of sight passing through the sphere. M is a suitable normalization 
constant. We have assumed here that single spheres are too small to be individually resolved. 



£ c and £ a are computed from Eq. (27) and Eq. (28), with E being the thermal energy. We 
find (Tb) = +1.22 niK at z = 25, increasing to a peak value of +1.41 mK at z = 22.7, and 
decreasing thereafter. The low value of (Tb) is due to the high degree of ionization of the 
central core and the rapid decrease in temperature away from the core. The temperature 
observed by an instrument = / (Tb) + (1 — /) T ot her- / is a "filling fraction" which takes 
into account that the individual spheres are in general, non-overlapping. T ot her represents 
the contribution due to other sources. Using the star formation rate of [46], we estimate 
/ ~ nR\n^(z = 25) ~ 0.3 where n*(z) is the number density of stars at z. Thus, if the dark 
matter is light (m x ~ 10 GeV), we expect gas heating from early star formation to be only 
a small correction to the heating caused by particle annihilation, at high redshifts. Heavier 
dark matter particles (m x ~ 50 GeV) with low concentration parameters C200 ~ 5 may be 
harder to distinguish from early stars. 



B. Global evolution of the brightness temperature. 

The hot bubble framework discussed in the previous subsection is not applicable at late 
times, when the bubbles coalesce. Let us briefly discuss the astrophysical models considered 
by Pritchard and Loeb [47] who compute the global evolution of the 21 cm brightness 
temperature in the presence of astrophysical sources. Three models (named A,B, and C) 
are considered in [47] depending on the star formation efficiency, and the number of ionizing, 
Ly-a, and X-ray photons produced by the sources. At a redshift z = 25, the spin temperature 
is nearly equal to the CMB temperature and the brightness temperature relative to the CMB 
~ -2 mK for all models (see Fig. 1 of [47]). For smaller redshifts, the spin temperature starts 
to follow the gas temperature again due to increased Ly-a coupling, resulting in a minimum 
at z ~ 17 with a H21 cm differential brightness temperature ~ —100 mK for models A 
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and B and ~ -75 mK for model C. The ionized fraction is still very small. The differential 
brightness temperature passes through zero, becoming positive at z ~ f 4 — f 5, as the gas is 
heated above the CMB temperature. A local maximum is reached at z ~ 12 — 14, quickly 
falling off to zero at lower redshifts as the Universe reionizes. 

These models based on gas heating by luminous sources also exhibit a minimum in the 
amplitude of fluctuations at redshifts z ~ 14 — 15, when the gas temperature exceeds the 
CMB temperature (as long still small), as described in the previous section. Thus 

detecting a minimum in the power spectrum at redshifts z < 15 cannot be interpreted as due 
to the presence of WIMP dark matter. For higher redshifts z > 20, astrophysical sources 
have very little influence upon the global gas temperature and the average 21 cm differential 
brightness temperature. We therefore expect our results based on dark matter heating to 
remain valid at these redshifts. The existence of a minimum in the amplitude of fluctuations 
at redshifts z > 20 would be a possible indication of heating by dark matter annihilation. 

VI. CONCLUSIONS. 

We studied the effect of dark matter annihilation on the CMB polarization and the H21 
cm power spectra, for different values of the particle mass, and assumed halo concentration 
parameter. We considered dark matter halos fitted with an NFW profile, and computed 
their luminosity. We used the Press-Schechter theory to determine the number of dark 
matter halos, and derived an expression for the energy absorbed by gas at a given redshift z. 
We compared this value to the energy generated at z and found that only a small fraction 
of the generated energy is absorbed (Fig. 1, peak absorption ~ 20% for 10 GeV WIMPs, 
and ~ 6% for 100 GeV WIMPs). About 30% of the absorbed energy goes into ionization 
and an equal amount into heating. We also considered the effect of varying the minimum 
halo mass which is determined by the particle physics model (Fig. 2). We did not include 
secondary ionizations, bremsstrahlung processes, etc, nor did we account for the increase in 
halo density with time. Thus, more accurate treatments could result in a larger absorption of 
energy. We computed the evolution of the ionized fraction (Fig. 3) and the gas temperature 
(Fig. 4) with redshift, for different values of the dark matter particle mass m x and assumed 
concentration parameter C2oo- 

We then investigated the observable effects of ionization. We assumed a reionization 
model characterized by a redshift such that x ion = 1 for z < z*. In models with dark 
matter annihilation, x ion ^ for z > z*. We compared these models with the standard reion- 
ization model and calculated the value of z* for different particle masses and concentration 
parameters, for an optical depth r = 0.087 (Table I). Figs. 5 and 6 show the reionization 
history and polarization power spectra for the model with m x = 1 GeV, C200 = 10. The 
change in the polarization power spectrum is small, but observable for very light (~ 1 GeV) 
dark matter particles, and for favorable halo parameters (C200 ~ 10). Future CMB experi- 
ments such as Planck and CMBPol may be able to distinguish the two scenarios, and hence 
distinguish light dark matter models from others. It will however be difficult to distinguish 
the dark matter reionization scenario from another model with gradual reionization. 

We also studied the effect of dark matter annihilation on the Hydrogen 21 cm power 
spectrum. In Fig. 7, we compared the collisional and Ly-a coupling terms. Fig. 8 shows the 
differential brightness temperature of H21 cm radiation T^, for m x = 10 and 50 GeV and 
for concentration parameters c 2 oo = 5, 10, and 20. Also shown is the case of no dark matter 
heating (C200 = 0). We then considered perturbations to the baryon density (S n ), ionized 
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fraction (8 X ), and gas temperature (5t), and computed the multifrequency angular power 
spectrum of fluctuations Ci(u, Au) in T^. Fig. 9 shows the time evolution and the relative 
amplitude of ionization and temperature fluctuations, assuming S x = 5^ = 5 n initially. It was 
shown that ionization and temperature fluctuations are not very significant at low redshifts 
(z < 30). We studied the variation of the power spectrum Ci(u,Au) with frequency Au. 
For I = 500, Ci » for Au > 1 MHz, while for I = 5000, Q falls off to w for Au > 0.1 
MHz. The frequency dependence of Ci(u,Au) may be used to separate the Hydrogen 21 
cm signal from the foreground contaminants. Fig. 11 shows the H21 power spectrum for 
different redshifts. Fig. 11(a) shows the case m x = 10 GeV, C200 — 5. In this model, the 
amplitude of fluctuations decreases as z is increased from 15 to 25, and then increases from 
25 to 30. This behavior is due to the gas temperature T approaching the CMB temperature 
T 7 and then increasing beyond T 7 . Fig. 11(b) (m x = 50 GeV, c 2 oo = 5) does not show a 
minimum at any redshift, and is a monotonous function of z. This is because for this model, 
the gas temperature T never exceeds the CMB temperature at any redshift in the considered 
range (compare Fig. 4(a) with Fig. 4(b)). Fig. 11(c) shows the power spectrum for the 
standard scenario with no dark matter heating. Fig. 12 shows the variation of the power 
spectrum with redshift in the range 15 to 30, for I = 1000, for different values of 0200- in 
(a), the power spectrum shows a minimum with redshift for all values of C200 7^ (compare 
with Fig. 11(a)). In (b), the models with C200 — 2.5 and 5 show no minimum (compare 
with Fig. 11(b)), while the model with C200 = 7.5 shows a minimum at z = 16.4. Future 
Hydrogen 21 cm experiments such as LOFAR and SKA may be able to distinguish light 
dark matter models (m x ~ 10 GeV) from the standard scenario by measuring the power 
spectrum at different redshifts. However, it would be difficult to distinguish models with 
larger dark matter masses (m x ~ 50 GeV) from the standard scenario, except for favorable 
halo parameters (C200 > 10). In certain dark matter models where the dark matter particle 
is produced by a non-thermal mechanism, the annihilation rate can be much larger than in 
the standard scenario considered here. In such cases, the ionization and heating rates are 
substantially higher [22] which may make it easier to distinguish the effects of dark matter 
annihilation from other sources. 

Finally, we discussed the effect of astrophysical sources on our calculations. We found 
that at high redshifts, heating by early stars is not a very large effect, especially for light 
dark matter particles. Later star formation at lower redshifts, and the formation of the first 
galaxies will likely be larger contaminants. We discussed the evolution of the brightness 
temperature at lower redshifts in the models of [47]. Heating by astrophysical sources 
also results in a minimum in the amplitude of fluctuations at redshifts z < 15. However 
for redshifts z > 20, we expect dark matter annihilation to be the relevant effect. Thus 
detecting a minimum in the amplitude of Hydrogen 21 cm fluctuations at z > 20 would be 
a possible indication of heating by dark matter annihilation. 
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